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AliSTRACT 


This paper gives a mathematically rigorous derivation of first-order 
corrections to multi- impulse approximations to tite solutions to space flight 
optimization problems with bang-bang control. The rocket is suoject to an 
inverse square gravitational force and to a thrust force with constant magni- 
tude. The mass decreases linearly with time. It is assumed that an optimal 
impulsive solution lias been obtained for a problem with given initial and final 
conditions. The method may then be used to obtain first-order corrections to 
the Initial values of the costate variables. Indications are given on how the 
theory may be extended to higher order corrections. The theory is applied to 
intercept and rendezvous problems. 


FIRST ORDLR IMPULSIVE SOLUTIONS 


J. F. Andrus* 

Univt?rsity of New Orleans, New Orleans, La. 

Introduction. This paper Is concerned primarily with first-order correc- 
tions to impulsive approximations to optimal space flight maneuvers. Theoretically, 
higher order corrections can be obtained in an analogous manner. An example 
problem is the determination of the ignition and burnout times and the tlme- 
histories of the thrust control angles which define a minimum time, exoatmos- 
pheric maneuver from a given initial orbit to a given final orbit. Several 
coast and powered phases may be Included. Only problems with conrtant thrust 
magnitude, F, and fuel burning rate magnitude, E, on each thrust arc are 
considered . 

The Impulsive solution is defined^ to be the limit of the bounded thrust 
solutions as 6 increases without bound, where F » cb and c is the constant 
exliaust velocity. The corrections are in terms of the parameter t ■ l/li. 

They are the first and higher order terms of the Taylor series expansions of 
the variables about e = 0. 


Research sponsored in part by Unlv. of New Orleans and by Northrop under 
NASA Contract NAS8-20082. 

Index categories: Spacecraft Navigation, Guidance, and FI Iglit-F’ath Control 
Systems. 

*Assoc. Prof, of Math. Formerly Member of Sr. Tech. Staff of Northrop in 
Huntsville, Ala. Member AL\A. 


1 


2 


The problem of obtaining corrections to Impulsive solutions has been 

2“6 

studied Independently by this author and several men associated with 
I’rinceton University. Refs (2) and (4) are concerned with the problem of 
constant F/m, whereas (3) and (5) extend the work to the problem considered 
in this paper. The latter works make use of expansions in terms of two 
parameters;namely, initial tlirust acceleration and the rocket Jet exliaust 
velocity. Ref. (6) considers applications to low-thrust mission analysis. 

The present paper is an extension of Ref. (7) to cover general initial 
conditions as well as final conditions. Moreover, the necessary conditions 
are given in a more explicit form. Both papers empiiasize mathematical rigor 
in the development of the corrections. 

State and Costate hquations. The basic problem is the determination of 
the optimal vacu im flight of a space vehicle from a prescribed set of Initial 
conditions co a given set of final conditions. The flight has at least one 
coasting piiase and is to be optimized with respect to payload. 

The equations of motion are 

y ” “ m - -B 

where k(A.) is the optimal steering vector \/|x|. The parameters F and B are 
zero on coast arcs. Tlie costate vector A is tlie solution to the costate 
equations 


• QU.y.A ). 


k 


Let tj^ and tj^ be the initial and 
1,2,...,N. Let < be the so-called 


final times on the k-th thrust arc for 
"switcliing function" satisfying the 


KKPHDDlU ll.l I. I I 


equation 
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- - U(A ) 
in 

where U - ^ A / 1 |x | . Tlie necessary conditions of optimality include the condi- 
tion ► ■ 0 at times tj^ and tj^ for k ■ 1,2,...,N except sometimes for time t^ 
and time tj^ if tliey corresj^ond to the initial time t^ and the final time t^, 
respectively. 

Impulsive Solutions. Assume that for each non-zero value of t In some 
neigliborhood of zero there is a solution y(t, ), y(t,t), A(t,f), X(t,t), ic(t,e), 
ra(t,t), t|^( ), tj^(i-), tj(t) to the boundary condition problem. ..'ere y, for 
example^ is cons’dered to be a function of the two arguments, t and e. It is 

also assumed that the impulsive solution exists; i.e., that the variables 
approach finite limits us l approaches zero. In Ref. (1) and later in this 
paper it is shown tliat, in the limit, t^^ “ t^^, y(tj^,0) “ y(tj^,0), A(tj^,0) « A(t|^,0), 

A(7j^.n) - A(tj^.O), 

_ m(t ,0) 

y(t ,0) = y(t. ,0) + c( n — ;; ) ,0)) 

m(tj.,0) 

m(t. ,0) 

'.V (0) = c -n 

m(tj^,0) 

where " y(t,^.0) i . 

The multi-point boundary-condition problem for tlie impulsive case involves 

the choice of y , y ,A ,A ,k , t , , t„ , . . . , t„, and t^ such that tc(t, ,0) = 0 for 
o ooooiz N 1 K 

k = 1,2,...,N; sucii that kU, .0) = 0 for k = 1,2,...,N (except sometimes for 

K 

k = 1 and/or k = N); and such that the given initial and final boundary conditions, 
including transversal 1 ty conditions and a scaling condition upon A, are all 


sat isf led . 


4 


Notation. Let yj^(t) - yltj^(0,c), “ y(t,^(e).£). y^ - 3y/3t, 


^ek 'kc^k* 


y^ “ y - 3y/3t, yj^^ 

Also let L*(t,e) ■ L[A(t,e)), C*(t,e) ■ G( t (c) , y ( t , e) ) , and so on. 

Tlie symbol o(c*^) will denote a finite summation of terms of the form 


n. n„ ^ n- 

a( )e "[At^(t)l ^(U^(e)l 


where nj^,n 2 ,n^, and n are non-negative integers such that nj^>n 2 +n^ > n and 
a(e) is bounded .;ithin some neighborhood of c ■ 0. Observe that U*(0) ■ 0 
whenever “ 0* The symbol is used Instead of when each 

coefficient a(c) in the summation is a constant. 

The General Procedure. Suppose for example that one of l 1 e constraints 
is y^ ■ a. Then y^(t) a is an identity in t. Therefore, yfj.(t) - Ot y^ ^ 

etc. If y^, corresponds to y^^, then 0* Tliis paper gives an expression 

which relates y, (0) to y. (0). Moreover, y. ™^y written in terms of 

KE iCC K » X I C 

yi, • .. and t. , . Ultimately the condition y^ (0) “ 0 can be expressed 

K f. K t K 1 K • X * C t C 

« ^ 

as a linear equation in the unknowns \ , \ , t, (0), etc. Similarly the other 

boundary conditions lead to linear equations. Tlie linear equations may be solved 
for the unknown derivatives. Once example lias been calculated, 

one may correct the Impulsive value ^^(0) hy adding the flrst-oruar correction 

EA (0) with £ >■ 1/l.A. 

OE 

Problems arise in the applications of the procedure Just described. These 
are elaborated upon below. 

Let iltj^ ° ^k ~ ^k ‘""k ™k ~ ""k’ “ ”^^k^ ’ L' Hospital's 

rule gives tlie equation Atj^ (0) » -Amj^(O). In general, 

I'l' 

. _j_^_ ^ 

, n n+1 , n+1 

ci> d ^ 
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at “0. Tlius, for example, once tj^ (C) 1 h known, t^^ (0) can be easily 
calculated. Kather tlian considering (0) to be the unknowns on 

the k-th thrust arc, the unknowns will be tj^^ (or tj^^) anc 

It will become evident that on any thrust arc for which *^1^(0) “ 0i 
conditions (0) ■ 0 and (0) ■ 0 are dependent. In order to obtain rn 

Independent condition, the condition “ 0 will be expanded and it will 

be discovered that the condition will involve only first dciivatives with respect 
to c. In general, the condition d" ■ 0 (or d"tCj^/dt" • 0) will be replaced 

u I'l'*"! ~ / j 0+1 _ 

by d 

Derivatives over Coast Arcs. In general 


yk+l.C ■ * l’>'k+l'®>'k“>'kc-'kc>'k> * l’>'k+l'*>’k'‘>'kc-'kES> 


■ 'kH..Cl * "\+l^’>'k'<>'kC-‘kE>'k> * ^’^tl'">'k'‘>'kE-'■k^‘^^ 


\+l,e ■ ‘k+l,c\+l ■" •“k+l'“>'k'<\c"'kt>'k> * '“k+l'’''k'<>'kc-‘kA> 


* I='*k+i^’*ki<‘kc-'kc‘k> * i’\*i'“'ki<’k -'kA» 


'kH.c • 'k+i.cVi * k 


k laitki'^\l<‘kc-'kk‘k> * l»kki'^»k)<>'k,-'kcV 


'\+l,c "'kc 

where the transition matrices of partial derivatives can often be expressed in 

2 

closed form as on a Kepler arc or In the case of G = y where .. is a const.iiu . 
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K 


ke 


The derivatives of ►: can be derived from the relationship “ 

+ Thus 


k+l,c “ “kc ^‘^^'"k^^‘k+l,L^k+l ■ *^kt^k 


* U/l\l)l>Ll“kU,c-'kU..Vu> - *I<V-'kc*k>>'- 


If “k ‘ 0 and i 0, 


the latter equation reduces to the simple condition 


^k+l^k+l,e‘ 


■ » 


( 1 ) 


applying wh 3n e ■ 0. 

Integrals which Approach Zero. Consider a function f (X ,X ,y ,y ,m, t) which 
has continuous partial derivatives and such that f*(t,e) ■ f [ A ( t , c ) , X ( t ,c ) ,y (t , e ) , 
y ( t , E ) ,m(t , c ) , t ] is continuous and bounded in tlie interval [tj^(E), tj^(c)] within 
some neighborhood of c ■ 0. Let 

T. (e) 

1(e) - -t (t) f*(t,e)dt 

By a mean value theorem for integrals, 

1(e) - f*[t(t),E) . dt - At (e)f*[t(E),cJ 

tk(e) k 

For some t in the interval [tj^.tj^]. Since f* is bounded, 1(e) » o(e). 

Now let 

.T 

■'k 


J(e) 


f*(t,E)dt dT 
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By the mean value theorem for double Integrula, 

J(r.) - f*(t,e) . dt dT ■ ■l/it^f * (t ,e ) - o(c^) 

"■k k 

for some t In the interval 

From the rules for differentiations under an Integral sign it follows 

that 

^ 1(0 ■ r, * /'l- tj dt 

k 


where 


■ f* dt 
'k ^ 




*T 

, (f, A. 

t, At 

k 


*T • 


+ f*^y + f*^y + f*m )dt 
e y c y t me 


If it is assumed that the components of f^ and the other partial derivatives 
are bounded in some neighborhood of c ■ 0, then the mean value theorem implies 
that 



* 

f dt 


c 




At. 


+ f (tm )] 
ir. t 


t-t 


where t is in [tj^.tj^). It must now be argued that the functions eA^[t(c),el. 

etc. are bounded for c 0 within some neighborhood of t • 0. For example, 

2 

it will be shown that cm ij bounded. On a thrust arc m • 1/c . Therefore 

e £ 


- t( )i l-(t)-tj^(L) 

tm^(t(£),£] - tm^(tj^(£),f.) -2 dt « rra^, ( tj^fc) , - ) f ^ 

k' c 


But ( t (• ) ]/t- J. tj^( ) where 11m l..tj_(>)/-] “ 

'0 


. t|^ (0). Therefore 
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in bounded within some nei({hburhuod of • ■ U. Also is bounded provided 

(in (t|^ bounded since un^ Is constant over a coast arc. Eventually the 

bounded’.tcss of (in (t, j depends upon whether or not (m (t , ) Is bounded; but 
c k ^ o 

in (t ,l) Is zero. Hence dl/di ■ j( ). 
o 


Similarly it can be argued that e'^3"m/3£'', e^S^y/Sc'^, etc. are bounded 
within some neighborhood of c - 0 for n - 1,2,... . The gr>eral Idea employed 
Is that factors of e In the denominator of the terms In the equations of variation 
are removed by multiplying by factors of c. Finally it can be concluded that 
£ d I/dE ■ o(£) and e d J/de ■ o(e). 


Now It Is possiole to define the symbol 
fl(c™), a functir>n ofe” ^)t'^d''l/de'^ (n ■ 0,1, 
(n ■ 0,1,2...), or a sum of such fur^tlonn. 
preceding paragraph It can be shown that 


Wfc™), m ^ 1, denoting a function 

\ £: ^ j / m™! V n*l .n , > I n 

...)^ a function o(£ )t d J/dc 

Employing the conclusions of the 


^ 0(£"’) - «(£"‘-P^ 
d£P 

for m > p > 0 and that 0((.) ■ o(e). Observe that one can not state with assurance 
that ^ o (c*") ■ oft'" V> because o(e) has a factor afc) which is bounded within 
some neighborhood of e ■ 0 but whose derivatives may not exist. Little is 
known about a'.e) because it may correspond to a function evaluated at an in- 
determinate time t(e). 

Changes over Thrust Arcs. In this section expressions of the form 
^ ^ ^ ^ ^ 0( ^) , . . . will be derived 

for the non- impulsive case. For example, the change in y over a thrust arc 
will be expressed as some function Yj^( ) plus a second-order term H((^). Similar 
expressions will be obtained for the first derivatives with respect to * . The 
exprecslons for tlie higher derivatives, though complicated, can be obtained 
similarly. 

lM:PK')l'l.d li.i ,1 i i U. . 

. IS 
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The first step Is to fir expressions for the Integrals of C* . Since 
•1/ , repeated integration by parts gives 

r t 

G*dt ■ ; mG* dt 

, 1 2, 2 *,^k 1 2, 2 * ,*/k . 3. 

- (mG*J - rp |m G I - T<^> Im G L 1 + «( ) 

^k ^ 'k ^ *^k 


and so on. Since the first derivatives will be of primary concern, ordinarily 

2 

the expansion will only be carried out to the second-order term Q(* ). The 

double integral is simply _ 

^ G*dt di - 4^ ^ lm^(J*)^*^ + eAt.m.Gf + «(c^) 

^ ^ 2 tj^ k k k 

Now the changes in y and y over a thrusting phase will be examined. Since 

y ■ (cr/m)L + G, 


c 1 * ^k I 

y “ fi y, ~ dt dl + .’ G*dt dx 

^k k^k . t^ m t^ tj^ 

^1. 

- - ctmj^fitn^ )Lj^ - cc(ml. ^ G* dt dt 

mj^ k ‘k *'k 

• c '"k 1 * ^k 

Ay, - - — dt + .• G* dt 

c t^ m tj^ 

- cfiln;^ )Lj^ + ccm^fJn — )L^ + cG(mI.j^]^ + ^ G* dt + 0(c ) 

m, m, k k k 

k k 


Since A - i^t.y.X), 


• L T 

A>, - At, A, + (}*dt dt 

k k k t, t, 

k k 




• *^k 

LA, - . Q* dt 

" ‘k 
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= - ELmO*].*" + • 

k 


Since k = - U(X ,X), 
m 


Ak, = c / ^ — dt = cell, £n — + ce^U m, £n — + Q(e^). 
k t, m k — t, K — 

k k 

An alternative procedure employing series expansions in terms of e leads 
to the same results. However, it is difficult or impossible to prove, in a 
logically rigorous way, that such series converge or that they can be 
differentiated. 

Changes over Thrust Arcs in the Limit. Letting e approach zero in the 
expressions obtained in the preceding equations, one obtains Ayj^(O) = 0, 

Ayk(O) = c £n(mj^/ d 1 j^)Lj^, AXj^(O) = 0, AXj^(O) = 0, and AKj^(O) = 0. Since |Lj^| = 1, 

it follows as in Ref. (1) that at e = 0, AVj^ = c £n(m^/ nij^) , 

— . * 

= i\ exp(-AVj^/c), and Ay^^ = AVj^Lj^. 

Taking derivatives of the expressions obtained for Ay^^, etc. in the 

preceding section and taking the limit as e apprc aches zero, one obtains 

— * 


Ay^^(O) - 


”k\ 

“k.«» - -"Vk- ' V\ 


2c „* 


''<kaa«> ■ — “k ^ “Vkc + 

W 


where = AVj_mj^ + Am^^c and a^^ = AV^^m^^ + Am^c. Since the direct calculation 
of aj^ and a^^ involves subtraction of nearly equal numbers, one should employ 
the series expansions 


(AV /c)^ Z 

k (n+2) : 
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Z (.V,/c)" 

n*0 

* 

If K - 0, then •• 0 so that 0* condition “ 0 

implies ■ 0. In this case tlie condition ■ 0 'ihould be emploved rather 

than one of the conditions, ^ and “ Ot in obtaining the first- 

order corrections. Tne condition ■ 0 is simply 


. X T* 

X, X, r X, X, 
k kc k lo: 


ih ^Vk ^ 

k 




( 2 ) 


Summary ot Conditions. The multi- point bounaary-conditlon problem will be 
summarized as an example for problems in which the initial and final phases of 
the flight are thrusting phases. 

T 

Corresponding to the initial point there is the scaling condition “ 1 

% 

and initial conditions h(yj^,yj^ .X^^.Xj^) “ 0, where h stands for a column vector 
with five or six components, some of which may represent transversallty con- 
ditions. (The conditions may be simply expressions of initial values of y and y.) 
The vector h has only ‘‘Ive components if tCj » 0, in which case condition (2) 
also applies. As e/.,.iained earlier the condition 0 is not employed. In 

any case the equatic'’ “ 0 need not be included explicitly because it can be 
trivially satisfied by setting “ 0. Therefore, corresponding to the initial 

point (or initial thrust arc), we nave the following seven conditions: 


T 


iih 


3h 


dh 


31) 


3yj^ ^Ic ^ 3yj^ ^It ^ 3X^ ^Ic 3X^^ ^le 




(applies if “ 0) 


J 

i{n’!Mii>rni:ii,rrv (jk uik 

■ ' J:’ 


V. 


There are also the 13 unkniwns, ^le’ ‘"it’ 

the Initial point. 

For each Intermediate thrust, conditions (1) and (2) apply and there 
'.re two unknivrns; namely, t^^ and on the k-th thrust arc. 

Condition (1) with k ■ N-1, corresponding to the equation ■ 0 Is 
T T — 

X A., - A^, A^, , “ 0. This corresponds to the final point. 

N Nt N-1 N-l,e ^ 

Corresponding to the final point chert* .are given conditions 

K(t.,,y.,,y.,,A.,, A^. k^,) * 0 where g has six or seven components, some of which 
N N N N N N 

may represent transversal Ity conditions. Tlie vector g has only six components 

In the case that ■ 0. Then the condition (2) with k » N applies. It 

corresponds to the condition - 0. In this case the condition k„ » 0 

Nee 2c 

or any equivalent condition must not be employed. Hence, corresponding to the 
final thrust arc, we have the eight conditions 


T T — 

^N^Ne ^N-l^N-l,e 


'T T* 

a‘a^ + 

N Ne N Nc 


4 ?; 


(applies i f • - 0) 


ju — 3t!— 3 b — Sc — *T 3u 

y , + ' -av„a.,a„ 

Nc — •^Mc — Ne — Ne N N — 

’S'n “n 


/ 


since “AV^A A . Tlie two unknowns, t„ and Am., , correspond to the final point, 
Nc N N N Nc Ne 

(If one of the conditions Is simply t„ = given numoer, then t =0.) 

N Ne 

Therefore, for the thrust- coast- thrust problem there ^re 15 linear equa- 
tions and 15 unknewns. For each additional intermediate thrust arc: there are 


two additlone’l equations and unknowns. It is understood that v., , for exatrnlc, 

• Ne 

in the above equations will be expressed, by means of equations derived earlier 


in this paper, in terms of the indicated unknowns. 
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An Intercept I'roblem. As an example consider .j tiirust-coast intercept 

problem in whicli the vehicle moves from a given point in state space to a 

point with given position components. The time of intercept is also specified. 

In this problem t must be identified with t, . Since A ■ 1. 

o 111’ 




(3) 


Since V- is fixed, we have v- - 0 so that 
2 •^2c 


_ + ^y 
^ y,i + — y 


3y 3y 


2 - + 
cl 

1 


But y t “ y, -t, y, 
•'el •'le le-'l 


— + — 


. + Am^y^ • and y^j - 


t, y, “Ay, + Am.G* 
Ic"^! ■'le I 1 


- — Am, A +AV A +a (I-A A/^)' at e “ 0. 
— lei lie 1 111 


Therefore, 


3y2 ^y, 

3y^ ^y^ 


— Am, A, + AV,A, + a, (1-A a/)A 

— lei llel 111 

mi 


0 (A) 


Equations (3) and (4) may be written in matrix form as 


c ^^2 

'l 

m^ 3y^ 




\ 


Ic 



0 
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where 



We observe that in general 


3y2 3y2 . 

X - ^ A + A 

3yi 3yi 


In an intercept problem the transversal ity conditions imply that A^ “ 0. Therefore, 




\^3yi J 3yi 


Utilizing the latter equation and equation (5), it may be sliown that 


Am 




Ic 


T* 


'u • ^ - »,l 


at t = 0. Starting with the requirement that A^ = 0, it nuiy also be shown 


that 
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Let tj^ ■ 0, t 2 ■ 380, u “ .388 x 10^^ m^/sec^, c ■ 4100 m/sec, 6 ■ 22 

kg-sec/m, ■ .168920 x 10^ kg-sec^/m, - (.287253 x 10^ m, .590785 x 10^ 

m, .777376 x 10^ m) , - (-7326.35 m/sec, 3219.14 m/sec, -474.472 m/sec), 

T 7 

y 2 * (0, .655630 x 10 m, 0). A solution to the boundary condition problem 
gives - 86.00 sec, - (.6618, .3727, .6505), - (-.1725 x lo”^, 

-.1198 X 10“^, -.1596 X lO"^). 

The impulsive solution yields tj^ “ ■ 76.64 X^^ - (.6631 , .3557, 

.6587), X^^ - (-. 1726 x lo"^, -.1149 x lo"^, -.1617 x lO"^) . A first-order 

— T 

correction to the impulsive solution gives t^^ “ 84.04, X^^ - (.6623, .3726, 

.6503), X^^ - (-.1726 x lO"^, -.1198 x lO"^, -.1596 x 10‘^). The error in 
Xj^^, for example, lias been reduced from (.00121, 01701, 00823) to (.00041, 

.00004, .00013). 

A Rendezvous Problem. Consider a thrust-coast-thrust rendezvous problem 

in which the vehicle moves from a given point in state space to point with 

given position and velocity components. The time of rendezvous is specified. 

For simplicity we take C -g where g Is constant. Tlierefore, () » X » 0. We 

identify t^ with t^^ and t^ with t 2 * The boundary conditions are y 2 (t) - constant, 
— T 

y 2 (c) = constant, X^ (c)X^(e) = l,»C 2 (f:) “ 0* 

Following a procedure similar to that of the intercept problem, we obtain 


Ic 


^ (ai + a2>.^^X2) 


m (a^-a_) 

‘”2£ - <‘i *2 - - '«1C 




/ . :v- M. PiHlR 


16 


Ic 


A - I- 

i\T Ic '> ’ 


a.-a„ _ 

^ (A, - (^2 X^IA^] 


At 2 ( ^^2 ‘ 


+ (m^ - 2nij^ + m^)[l - (A^ Xj^lJX^ 


7 

ii 


where AT ■ t2 ~ tj* 

Conaider the numerical problem in which tj^ ■ 0, ^2 “ sec, " .17 

5 2 T 2 2 

X 10 kg-sec /m, ti - 22 kg-sec/m, c * 4100 m/sec, g - (4.63 m/sec , 8.00 m/sec ), 

- (.18 X 10^ m, .63 x 10^ m) , y 2 ^ ■ (.413102 x 10^ m, .501462 x 10^ m) , 

y^^ - (.68 X 10^ m/sec, -.2 x 10^ m/sec), y 2 ^ * (.528352 x 10^ m/sec, -.505836 

4 • T 

X 10 m/sec). The optimum solution is tj - 50 sec, t 2 “ 350, A^ •• (.8, .6), 

aJ - (-.002, -.004). 

The impulsive solution yields t^^ = 46.47, t 2 = 352.47, = (.8414, .5404), 

Aj^ « (-.001733, -.004010). Tlie impulsive solution witli a first order correction 
added is 7^ - 49.64, t 2 - 350.60, A^^ » (.7992, .6060), kj - (-.001979, -.004063). 
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